act_site = [1 2 4 7 8];
pas_site = setdiff(1:9,act_site);
con_site = 10;

for time = 1:7
    for site = 1:9
        cur_vel = -proj_meta(site).rd(1,time).velM_smoothed(:,7501:15000);
        cur_vel(cur_vel > 0.005) = 1;
        cur_vel(cur_vel <= 0.005) = 0;
        run_time{site,time} = sum(cur_vel)/length(cur_vel)*100;
    end
    run_time_act_m(time) = mean([run_time{act_site,time}]);
    run_time_pas_m(time) = mean([run_time{pas_site,time}]);
    run_time_act_s(time) = std([run_time{act_site,time}])/sqrt(5);
    run_time_pas_s(time) = std([run_time{pas_site,time}])/sqrt(4);
end

figure;hold on
errorbar(run_time_act_m,run_time_act_s,'--')
errorbar(run_time_pas_m,run_time_pas_s,'r--')